Integrating multiplexed imaging and multiscale modeling identifies tumor phenotype conversion as a critical component of therapeutic T cell efficacy.

Cancer progression is a complex process involving interactions that unfold across molecular, cellular, and tissue scales. These multiscale interactions have been difficult to measure and to simulate. Here, we integrated CODEX multiplexed tissue imaging with multiscale modeling software to model key action points that influence the outcome of T cell therapies with cancer. The initial phenotype of therapeutic T cells influences the ability of T cells to convert tumor cells to an inflammatory, anti-proliferative phenotype. This T cell phenotype could be preserved by structural reprogramming to facilitate continual tumor phenotype conversion and killing. One takeaway is that controlling the rate of cancer phenotype conversion is critical for control of tumor growth. The results suggest new design criteria and patient selection metrics for T cell therapies, call for a rethinking of T cell therapeutic implementation, and provide a foundation for synergistically integrating multiplexed imaging data with multiscale modeling of the cancer-immune interface. A record of this paper's transparent peer review process is included in the supplemental information.

Spatially Segregated Macrophage Populations Predict Distinct Outcomes In Colon Cancer.

Tumor-associated macrophages are transcriptionally heterogeneous, but the spatial distribution and cell interactions that shape macrophage tissue roles remain poorly characterized. Here, we spatially resolve five distinct human macrophage populations in normal and malignant human breast and colon tissue and reveal their cellular associations. This spatial map reveals that distinct macrophage populations reside in spatially segregated micro-environmental niches with conserved cellular compositions that are repeated across healthy and diseased tissue. We show that IL4I1+ macrophages phagocytose dying cells in areas with high cell turnover and predict good outcome in colon cancer. In contrast, SPP1+ macrophages are enriched in hypoxic and necrotic tumor regions and portend worse outcome in colon cancer. A subset of FOLR2+ macrophages is embedded in plasma cell niches. NLRP3+ macrophages co-localize with neutrophils and activate an inflammasome in tumors. Our findings indicate that a limited number of unique human macrophage niches function as fundamental building blocks in tissue.

Subcortical-cortical white matter connectivity in adults with autism spectrum disorder and schizophrenia patients.

Autism spectrum disorder (ASD) and schizophrenia (SZ) are neuropsychiatric disorders that overlap in symptoms associated with social-cognitive impairment. Alterations of the cingulate cortex, subcortical, medial-temporal, and orbitofrontal structures are frequently reported in both disorders. In this study, we examined white-matter connectivity between these structures in adults with ASD and SZ patients compared with their respective neurotypical controls and indirectly with each other, using probabilistic and local DTI tractography. This exploratory study utilized publicly available neuroimaging databases, of adults with ASD (ABIDE II; n = 28) and SZ (COBRE; n = 38), age-gender matched neurotypicals (NT) and associated phenotypic data. Tractography was performed using Freesurfer and MRtrix software, and diffusion metrics of white-matter tracts between cingulate-, orbitofrontal- cortices, subcortical structures, parahippocampal, entorhinal cortex were assessed. In ASD, atypical diffusivity parameters were found in the isthmus cingulate and parahippocampal connectivity to subcortical and rostral-anterior cingulate, which were also associated with IQ and social skills (SRS). In contrast, atypical diffusivity parameters were observed between the medial-orbitofrontal cortex and subcortical structures in SZ, and were associated with executive function (i.e., IQ, processing speed) and emotional regulation. Overall, the results suggest that defects in the isthmus cingulate, medial-orbitofrontal, and striato-limbic white matter connectivity may help unravel the neural underpinnings of executive and social-emotional dysfunction at the core of neuropsychiatric disorders.

Epstein-Barr Virus Orchestrates Spatial Reorganization and Immunomodulation within the Classic Hodgkin Lymphoma Tumor Microenvironment.

Classic Hodgkin Lymphoma (cHL) is a tumor composed of rare malignant Hodgkin and Reed-Sternberg (HRS) cells nested within a T-cell rich inflammatory immune infiltrate. cHL is associated with Epstein-Barr Virus (EBV) in 25% of cases. The specific contributions of EBV to the pathogenesis of cHL remain largely unknown, in part due to technical barriers in dissecting the tumor microenvironment (TME) in high detail. Herein, we applied multiplexed ion beam imaging (MIBI) spatial pro-teomics on 6 EBV-positive and 14 EBV-negative cHL samples. We identify key TME features that distinguish between EBV-positive and EBV-negative cHL, including the relative predominance of memory CD8 T cells and increased T-cell dysfunction as a function of spatial proximity to HRS cells. Building upon a larger multi-institutional cohort of 22 EBV-positive and 24 EBV-negative cHL samples, we orthogonally validated our findings through a spatial multi-omics approach, coupling whole transcriptome capture with antibody-defined cell types for tu-mor and T-cell populations within the cHL TME. We delineate contrasting transcriptomic immunological signatures between EBV-positive and EBV-negative cases that differently impact HRS cell proliferation, tumor-immune interactions, and mecha-nisms of T-cell dysregulation and dysfunction. Our multi-modal framework enabled a comprehensive dissection of EBV-linked reorganization and immune evasion within the cHL TME, and highlighted the need to elucidate the cellular and molecular fac-tors of virus-associated tumors, with potential for targeted therapeutic strategies.

A Spatial Multi-Modal Dissection of Host-Microbiome Interactions within the Colitis Tissue Microenvironment.

The intricate and dynamic interactions between the host immune system and its microbiome constituents undergo dynamic shifts in response to perturbations to the intestinal tissue environment. Our ability to study these events on the systems level is significantly limited by in situ approaches capable of generating simultaneous insights from both host and microbial communities. Here, we introduce Microbiome Cartography (MicroCart), a framework for simultaneous in situ probing of host features and its microbiome across multiple spatial modalities. We demonstrate MicroCart by comprehensively investigating the alterations in both gut host and microbiome components in a murine model of colitis by coupling MicroCart with spatial proteomics, transcriptomics, and glycomics platforms. Our findings reveal a global but systematic transformation in tissue immune responses, encompassing tissue-level remodeling in response to host immune and epithelial cell state perturbations, and bacterial population shifts, localized inflammatory responses, and metabolic process alterations during colitis. MicroCart enables a deep investigation of the intricate interplay between the host tissue and its microbiome with spatial multiomics.

Multimodal and spatially resolved profiling identifies distinct patterns of T cell infiltration in nodal B cell lymphoma entities.

The redirection of T cells has emerged as an attractive therapeutic principle in B cell non-Hodgkin lymphoma (B-NHL). However, a detailed characterization of lymphoma-infiltrating T cells across B-NHL entities is missing. Here we present an in-depth T cell reference map of nodal B-NHL, based on cellular indexing of transcriptomes and epitopes, T cell receptor sequencing, flow cytometry and multiplexed immunofluorescence applied to 101 lymph nodes from patients with diffuse large B cell, mantle cell, follicular or marginal zone lymphoma, and from healthy controls. This multimodal resource revealed quantitative and spatial aberrations of the T cell microenvironment across and within B-NHL entities. Quantitative differences in PD1+ TCF7- cytotoxic T cells, T follicular helper cells or IKZF3+ regulatory T cells were linked to their clonal expansion. The abundance of PD1+ TCF7- cytotoxic T cells was associated with poor survival. Our study portrays lymphoma-infiltrating T cells with unprecedented comprehensiveness and provides a unique resource for the investigation of lymphoma biology and prognosis.

MAPS: pathologist-level cell type annotation from tissue images through machine learning.

Highly multiplexed protein imaging is emerging as a potent technique for analyzing protein distribution within cells and tissues in their native context. However, existing cell annotation methods utilizing high-plex spatial proteomics data are resource intensive and necessitate iterative expert input, thereby constraining their scalability and practicality for extensive datasets. We introduce MAPS (Machine learning for Analysis of Proteomics in Spatial biology), a machine learning approach facilitating rapid and precise cell type identification with human-level accuracy from spatial proteomics data. Validated on multiple in-house and publicly available MIBI and CODEX datasets, MAPS outperforms current annotation techniques in terms of speed and accuracy, achieving pathologist-level precision even for typically challenging cell types, including tumor cells of immune origin. By democratizing rapidly deployable and scalable machine learning annotation, MAPS holds significant potential to expedite advances in tissue biology and disease comprehension.

BCG vaccination stimulates integrated organ immunity by feedback of the adaptive immune response to imprint prolonged innate antiviral resistance.

Bacille Calmette-Guérin (BCG) vaccination can confer nonspecific protection against heterologous pathogens. However, the underlying mechanisms remain mysterious. We show that mice vaccinated intravenously with BCG exhibited reduced weight loss and/or improved viral clearance when challenged with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2 B.1.351) or PR8 influenza. Protection was first evident between 14 and 21 d post-vaccination and lasted ∼3 months. Notably, BCG induced a biphasic innate response and robust antigen-specific type 1 helper T cell (TH1 cell) responses in the lungs. MyD88 signaling was essential for innate and TH1 cell responses, and protection against SARS-CoV-2. Depletion of CD4+ T cells or interferon (IFN)-γ activity before infection obliterated innate activation and protection. Single-cell and spatial transcriptomics revealed CD4-dependent expression of IFN-stimulated genes in lung myeloid and epithelial cells. Notably, BCG also induced protection against weight loss after mouse-adapted SARS-CoV-2 BA.5, SARS-CoV and SHC014 coronavirus infections. Thus, BCG elicits integrated organ immunity, where CD4+ T cells feed back on tissue myeloid and epithelial cells to imprint prolonged and broad innate antiviral resistance.

T cell-mediated curation and restructuring of tumor tissue coordinates an effective immune response.

Antigen-specific T cells traffic to, are influenced by, and create unique cellular microenvironments. Here we characterize these microenvironments over time with multiplexed imaging in a melanoma model of adoptive T cell therapy and human patients with melanoma treated with checkpoint inhibitor therapy. Multicellular neighborhood analysis reveals dynamic immune cell infiltration and inflamed tumor cell neighborhoods associated with CD8+ T cells. T cell-focused analysis indicates T cells are found along a continuum of neighborhoods that reflect the progressive steps coordinating the anti-tumor immune response. More effective anti-tumor immune responses are characterized by inflamed tumor-T cell neighborhoods, flanked by dense immune infiltration neighborhoods. Conversely, ineffective T cell therapies express anti-inflammatory cytokines, resulting in regulatory neighborhoods, spatially disrupting productive T cell-immune and -tumor interactions. Our study provides in situ mechanistic insights into temporal tumor microenvironment changes, cell interactions critical for response, and spatial correlates of immunotherapy outcomes, informing cellular therapy evaluation and engineering.

Integrating Multiplexed Imaging and Multiscale Modeling Identifies Tumor Phenotype Transformation as a Critical Component of Therapeutic T Cell Efficacy.

Cancer progression is a complex process involving interactions that unfold across molecular, cellular, and tissue scales. These multiscale interactions have been difficult to measure and to simulate. Here we integrated CODEX multiplexed tissue imaging with multiscale modeling software, to model key action points that influence the outcome of T cell therapies with cancer. The initial phenotype of therapeutic T cells influences the ability of T cells to convert tumor cells to an inflammatory, anti-proliferative phenotype. This T cell phenotype could be preserved by structural reprogramming to facilitate continual tumor phenotype conversion and killing. One takeaway is that controlling the rate of cancer phenotype conversion is critical for control of tumor growth. The results suggest new design criteria and patient selection metrics for T cell therapies, call for a rethinking of T cell therapeutic implementation, and provide a foundation for synergistically integrating multiplexed imaging data with multiscale modeling of the cancer-immune interface.

Integration of spatial and single-cell data across modalities with weakly linked features.

Although single-cell and spatial sequencing methods enable simultaneous measurement of more than one biological modality, no technology can capture all modalities within the same cell. For current data integration methods, the feasibility of cross-modal integration relies on the existence of highly correlated, a priori 'linked' features. We describe matching X-modality via fuzzy smoothed embedding (MaxFuse), a cross-modal data integration method that, through iterative coembedding, data smoothing and cell matching, uses all information in each modality to obtain high-quality integration even when features are weakly linked. MaxFuse is modality-agnostic and demonstrates high robustness and accuracy in the weak linkage scenario, achieving 20~70% relative improvement over existing methods under key evaluation metrics on benchmarking datasets. A prototypical example of weak linkage is the integration of spatial proteomic data with single-cell sequencing data. On two example analyses of this type, MaxFuse enabled the spatial consolidation of proteomic, transcriptomic and epigenomic information at single-cell resolution on the same tissue section.

MAPS: Pathologist-level cell type annotation from tissue images through machine learning.

Highly multiplexed protein imaging is emerging as a potent technique for analyzing protein distribution within cells and tissues in their native context. However, existing cell annotation methods utilizing high-plex spatial proteomics data are resource intensive and necessitate iterative expert input, thereby constraining their scalability and practicality for extensive datasets. We introduce MAPS (Machine learning for Analysis of Proteomics in Spatial biology), a machine learning approach facilitating rapid and precise cell type identification with human-level accuracy from spatial proteomics data. Validated on multiple in-house and publicly available MIBI and CODEX datasets, MAPS outperforms current annotation techniques in terms of speed and accuracy, achieving pathologist-level precision even for challenging cell types, including tumor cells of immune origin. By democratizing rapidly deployable and scalable machine learning annotation, MAPS holds significant potential to expedite advances in tissue biology and disease comprehension.

Expanded vacuum-stable gels for multiplexed high-resolution spatial histopathology.

Cellular organization and functions encompass multiple scales in vivo. Emerging high-plex imaging technologies are limited in resolving subcellular biomolecular features. Expansion Microscopy (ExM) and related techniques physically expand samples for enhanced spatial resolution, but are challenging to be combined with high-plex imaging technologies to enable integrative multiscaled tissue biology insights. Here, we introduce Expand and comPRESS hydrOgels (ExPRESSO), an ExM framework that allows high-plex protein staining, physical expansion, and removal of water, while retaining the lateral tissue expansion. We demonstrate ExPRESSO imaging of archival clinical tissue samples on Multiplexed Ion Beam Imaging and Imaging Mass Cytometry platforms, with detection capabilities of > 40 markers. Application of ExPRESSO on archival human lymphoid and brain tissues resolved tissue architecture at the subcellular level, particularly that of the blood-brain barrier. ExPRESSO hence provides a platform for extending the analysis compatibility of hydrogel-expanded biospecimens to mass spectrometry, with minimal modifications to protocols and instrumentation.

Organization of the human intestine at single-cell resolution.

The intestine is a complex organ that promotes digestion, extracts nutrients, participates in immune surveillance, maintains critical symbiotic relationships with microbiota and affects overall health1. The intesting has a length of over nine metres, along which there are differences in structure and function2. The localization of individual cell types, cell type development trajectories and detailed cell transcriptional programs probably drive these differences in function. Here, to better understand these differences, we evaluated the organization of single cells using multiplexed imaging and single-nucleus RNA and open chromatin assays across eight different intestinal sites from nine donors. Through systematic analyses, we find cell compositions that differ substantially across regions of the intestine and demonstrate the complexity of epithelial subtypes, and find that the same cell types are organized into distinct neighbourhoods and communities, highlighting distinct immunological niches that are present in the intestine. We also map gene regulatory differences in these cells that are suggestive of a regulatory differentiation cascade, and associate intestinal disease heritability with specific cell types. These results describe the complexity of the cell composition, regulation and organization for this organ, and serve as an important reference map for understanding human biology and disease.

Organ Mapping Antibody Panels: a community resource for standardized multiplexed tissue imaging.

Multiplexed antibody-based imaging enables the detailed characterization of molecular and cellular organization in tissues. Advances in the field now allow high-parameter data collection (>60 targets); however, considerable expertise and capital are needed to construct the antibody panels employed by these methods. Organ mapping antibody panels are community-validated resources that save time and money, increase reproducibility, accelerate discovery and support the construction of a Human Reference Atlas.

Concerted epithelial and stromal changes during progression of Barrett's Esophagus to invasive adenocarcinoma exposed by multi-scale, multi-omics analysis.

Esophageal adenocarcinoma arises from Barrett's esophagus, a precancerous metaplastic replacement of squamous by columnar epithelium in response to chronic inflammation. Multi-omics profiling, integrating single-cell transcriptomics, extracellular matrix proteomics, tissue-mechanics and spatial proteomics of 64 samples from 12 patients' paths of progression from squamous epithelium through metaplasia, dysplasia to adenocarcinoma, revealed shared and patient-specific progression characteristics. The classic metaplastic replacement of epithelial cells was paralleled by metaplastic changes in stromal cells, ECM and tissue stiffness. Strikingly, this change in tissue state at metaplasia was already accompanied by appearance of fibroblasts with characteristics of carcinoma-associated fibroblasts and of an NK cell-associated immunosuppressive microenvironment. Thus, Barrett's esophagus progresses as a coordinated multi-component system, supporting treatment paradigms that go beyond targeting cancerous cells to incorporating stromal reprogramming.

tidytof: a user-friendly framework for scalable and reproducible high-dimensional cytometry data analysis.

Summary: While many algorithms for analyzing high-dimensional cytometry data have now been developed, the software implementations of these algorithms remain highly customized-this means that exploring a dataset requires users to learn unique, often poorly interoperable package syntaxes for each step of data processing. To solve this problem, we developed {tidytof}, an open-source R package for analyzing high-dimensional cytometry data using the increasingly popular 'tidy data' interface. Availability and implementation: {tidytof} is available at https://github.com/keyes-timothy/tidytof and is released under the MIT license. It is supported on Linux, MS Windows and MacOS. Additional documentation is available at the package website (https://keyes-timothy.github.io/tidytof/). Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Dasatinib overcomes glucocorticoid resistance in B-cell acute lymphoblastic leukemia.

Resistance to glucocorticoids (GC) is associated with an increased risk of relapse in B-cell progenitor acute lymphoblastic leukemia (BCP-ALL). Performing transcriptomic and single-cell proteomic studies in healthy B-cell progenitors, we herein identify coordination between the glucocorticoid receptor pathway with B-cell developmental pathways. Healthy pro-B cells most highly express the glucocorticoid receptor, and this developmental expression is conserved in primary BCP-ALL cells from patients at diagnosis and relapse. In-vitro and in vivo glucocorticoid treatment of primary BCP-ALL cells demonstrate that the interplay between B-cell development and the glucocorticoid pathways is crucial for GC resistance in leukemic cells. Gene set enrichment analysis in BCP-ALL cell lines surviving GC treatment show enrichment of B cell receptor signaling pathways. In addition, primary BCP-ALL cells surviving GC treatment in vitro and in vivo demonstrate a late pre-B cell phenotype with activation of PI3K/mTOR and CREB signaling. Dasatinib, a multi-kinase inhibitor, most effectively targets this active signaling in GC-resistant cells, and when combined with glucocorticoids, results in increased cell death in vitro and decreased leukemic burden and prolonged survival in an in vivo xenograft model. Targeting the active signaling through the addition of dasatinib may represent a therapeutic approach to overcome GC resistance in BCP-ALL.

Physioxia improves the selectivity of hematopoietic stem cell expansion cultures.

Hematopoietic stem cells (HSCs) are a rare type of hematopoietic cell that can entirely reconstitute the blood and immune system after transplantation. Allogeneic HSC transplantation (HSCT) is used clinically as a curative therapy for a range of hematolymphoid diseases; however, it remains a high-risk therapy because of its potential side effects, including poor graft function and graft-versus-host disease (GVHD). Ex vivo HSC expansion has been suggested as an approach to improve hematopoietic reconstitution in low-cell dose grafts. Here, we demonstrate that the selectivity of polyvinyl alcohol (PVA)-based mouse HSC cultures can be improved using physioxic culture conditions. Single-cell transcriptomic analysis helped confirm the inhibition of lineage-committed progenitor cells in physioxic cultures. Long-term physioxic expansion also afforded culture-based ex vivo HSC selection from whole bone marrow, spleen, and embryonic tissues. Furthermore, we provide evidence that HSC-selective ex vivo cultures deplete GVHD-causing T cells and that this approach can be combined with genotoxic-free antibody-based conditioning HSCT approaches. Our results offer a simple approach to improve PVA-based HSC cultures and the underlying molecular phenotype, and highlight the potential translational implications of selective HSC expansion systems for allogeneic HSCT.

A tissue atlas of ulcerative colitis revealing evidence of sex-dependent differences in disease-driving inflammatory cell types and resistance to TNF inhibitor therapy.

Although literature suggests that resistance to TNF inhibitor (TNFi) therapy in patients with ulcerative colitis (UC) is partially linked to immune cell populations in the inflamed region, there is still substantial uncertainty underlying the relevant spatial context. Here, we used the highly multiplexed immunofluorescence imaging technology CODEX to create a publicly browsable tissue atlas of inflammation in 42 tissue regions from 29 patients with UC and 5 healthy individuals. We analyzed 52 biomarkers on 1,710,973 spatially resolved single cells to determine cell types, cell-cell contacts, and cellular neighborhoods. We observed that cellular functional states are associated with cellular neighborhoods. We further observed that a subset of inflammatory cell types and cellular neighborhoods are present in patients with UC with TNFi treatment, potentially indicating resistant niches. Last, we explored applying convolutional neural networks (CNNs) to our dataset with respect to patient clinical variables. We note concerns and offer guidelines for reporting CNN-based predictions in similar datasets.